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High abundances of iodine monoxide (IO) are known to exist and to participate in local photochem- 
istry of the marine boundary layer. Of particular interest are the roles IO plays in the formation of 
new particles in coastal marine environments and in depletion episodes of ozone and mercury in 
the Arctic polar spring. This paper describes a ground-based instrument that measures IO at mix- 
ing ratios less than one part in 10 12 . The IO radical is measured by detecting laser-induced fluores- 
cence at wavelengths longer that 500 nm. Tunable visible light is used to pump the A 2 113/2 1 v' = 2) 
X 2 n 3/2 (v" = 0) transition of IO near 445 nm. The laser light is produced by a solid-state, 
Nd:YAG-pumped Ti:Sapphire laser at 5 kHz repetition rate. The laser-induced fluorescence instru- 
ment performs reliably with very high signal-to-noise ratios (> 10 ) achieved in short integration times 
(<1 min). The observations from a validation deployment to the Shoals Marine Lab on Appledore 
Island, ME are presented and are broadly consistent with in situ observations from European Coastal 
Sites. Mixing ratios ranged from the instrumental detection limit (<1 pptv) to 10 pptv. These data 
represent the first in situ point measurements of IO in North America. © 2014 AIP Publishing LLC. 
[http://dx.doi.org/10.1063/L4869857] 


I. INTRODUCTION 

Recent studies of atmospheric iodine suggest that it is 
important to broad ranging aspects of atmospheric chemistry 
with significant potential impacts on human health and cli- 
mate. Figure 1 highlights the various roles of iodine photo- 
chemistry in the atmosphere. First, laboratory studies 1-4 and 
atmospheric observations 5-13 imply that iodine contributes to 
the formation and growth of particles in the marine boundary 
layer (MBL). Second, atmospheric iodine may facilitate the 
oxidation of elemental mercury to a bioavailable form that can 
accumulate in ecosystems . 16, 17 Finally, reactive iodine can in- 
fluence the oxidative capacity of the atmosphere through the 
catalytic destruction of ozone 1 4-21 and the conversion of HCL 
to OH . 14 - 22 - 2 ’ 44 A recent review by Saiz-Lopez et al. 24 de- 
tails the current state of understanding of atmospheric iodine 
chemistry. 

The sources of reactive iodine in the atmosphere are 
short-lived iodocarbons (e.g., CH 3 I and CHiClI) and molecu- 
lar iodine, which have mixing ratios of a few parts per trillion 
(pptv, equivalent to pmol mol -1 ) in the MBL . 13,25 22 In recent 
years, emissions of reactive iodine precursors from macroal- 
gae in coastal regions have been studied extensively , 13, 26,29-31 
but these sources likely represent a small fraction of emis- 
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sions on the global scale . 24 In contrast, sources of organic and 
inorganic iodine over the open ocean are neither well under- 
stood nor well quantified . 24 Recent laboratory studies suggest 
that hypoidous acid (HOI) is an important source of atmo- 
spheric iodine , 32 but current analytical techniques are inade- 
quate for measuring ambient concentrations of HOI. Once in 
the atmosphere, these very short-lived sources of iodine are 
readily photolyzed to form atomic iodine, which then reacts 


with O 3 to form iodine monoxide (IO): 


RI + hv R + I, 

(la) 

I 2 + hv -* I + 1 , 

(lb) 

I + O 3 — ► IO + 0 2 ■ 

( 2 ) 


These reactions are known to produce significant amounts of 
reactive iodine (pptv of IO, OIO, HOI) in the MBL, " ,34, 41 
but longer-lived species (e.g., CH 3 I) may be transported to 
the higher altitudes by convective transport . 20,35 

Despite the potentially important roles played by 
IO in the various facets of atmospheric photochemistry, 
the database of available measurements to character- 
ize the distribution of reactive iodine is limited. The 
majority of measurements have been made using spa- 
tially averaged techniques based on differential optical 
absorption spectroscopy (DOAS), including long-path 
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FIG. 1 . The different facets of iodine photochemistry in the atmosphere ini- 
tiated from very short-lived source gases (VSLS): (A) ozone depletion, (B) 
Arctic Mercury Depletion Episodes (AMDE), and (C) coastal new particle 
and cloud condensation nuclei (CCN) formation. The transport of iodine 
across the tropopause (TTT = Thermal Tropical Tropopause, TTL = Trop- 
ical Tropopause Layer) does not lead to widespread elevation in IO concen- 
trations in the stratosphere. 


(LP-DOAS ) 6,9,13,14,33,34,36-40 and multi-axis (MAX-DOAS) 
observations . 42-47 Available in situ measurements by laser- 
induced fluorescence (LIF) 8 - U) - 4S and cavity ringdown 
spectroscopy (CRDS ) 7 show significantly higher maximal 
concentrations of IO than spatially averaged measurements. 
In an instrumental intercomparison, Commane et a /. 48 
reported that LIF measurements were as much as an order 
of magnitude higher than the DOAS measurements. The 
presence of “hot spots” that are highly episodic and spatially 
inhomogeneous suggests that in situ observations are critical 
for understanding the relevance of IO to nucleation events 1 
and radical chemistry, which may be sensitive to variations in 
abundance on short spatial scales. 

A global-scale climatology of reactive iodine requires 
a diverse database of measurements. Most direct observa- 
tions of reactive iodine and its precursors have occurred in 
a relatively narrow range of mid latitude coastal sites (e.g., 
Mace Head, Ireland and Roscoff, France). Ship-based ob- 
servations by MAX-DOAS have only recently been reported 
from the eastern 49 and western 43 Pacific Ocean. The satellite- 
based Scanning Imaging Absorption Spectrometer for Atmo- 
spheric Chartography (SCIAMACHY) provided global cov- 
erage, but absolute column densities and spatial distributions 
proved to be sensitive to the assumptions of the retrieval al- 
gorithms used to interpret the observations . 50,51 Direct, ver- 
tically resolved measurements of IO are necessary to pro- 
vide constraints for improving future retrievals of satellite 
observations. 


Highly sensitive in situ measurements are also critical for 
probing local concentrations of IO to determine its contribu- 
tion to ozone destruction throughout the free troposphere and 
lower stratosphere. The abundance of CH 3 I, a relatively stable 
photochemical precursor of IO, is known to be variable above 
the boundary layer and strongly dependent on local convec- 
tion (e.g., Refs. 52 and 53), which suggests that IO concen- 
trations may be highly variable as well. IO is not ubiquitous 
in the lower stratosphere 54-57 but may be significant locally . 49 
Though no in situ measurements in the free troposphere exist, 
MAX DOAS measurements of the slant column density of IO 
have been made recently at a high altitude observatory 47 and 
by aircraft. 4 "’ Both observations show detectable levels of IO 
at high altitude. Dix et a /. 45 suggest that more than half of the 
vertical column abundance of IO may exist above the bound- 
ary layer in the tropical free troposphere. The possibility that 
IO may be widespread throughout the free troposphere is sig- 
nificant: (1) It suggests that the sources and sinks of IO may 
not be fully understood; (2) It challenges the assumption that 
IO is localized to the boundary layer, which serves as a ba- 
sis for the interpretation of satellite observations ; 50,51 (3) It 
suggests that IO may influence tropospheric ozone. A mod- 
eling study based on observations of reactive halogens and 
precursors suggested the halogen-mediated ozone depletion 
accounts for 15%-30% of ozone loss integrated over the to- 
tal tropospheric column in the tropics . 39 Inclusion of iodine 
chemistry in the model led to a 4-fold increase in halogen- 
mediated ozone loss relative to bromine alone. These observa- 
tions and modeling results motivate future vertically resolved 
measurements of IO in situ throughout the free troposphere 
and into the lower stratosphere. 

The inhomogeneity of IO in the atmosphere requires a 
mobile in situ instrument to target regions that contain sig- 
nificant IO photochemistry. This paper describes the develop- 
ment and calibration of an IO instrument based on LIF de- 
tection, which has been used extensively by this laboratory 
for measuring free radicals in situ in the free troposphere and 
stratosphere . 98-61 The IO instrument was designed and tested 
for ground-based measurements with the capability of becom- 
ing deployable in future aircraft- and ship-borne campaigns. 
Preliminary measurements taken at the Shoals Marine Labo- 
ratory on Appledore Island, ME are presented. IO measure- 
ments with sub-pptv sensitivity are demonstrated with short 
integration times (~1 min). 

II. LIF AND THE DETECTION OF ATMOSPHERIC IO 

The detection of IO by LIF has been successfully demon- 
strated by several groups (e.g.. Refs. 8 , 62-65). The (v', 0), v' 
= 0, 1, 3, 4, and 5, bands of the A 2 n 3/2 <— X 2 n 3/2 transi- 
tion of IO do not exhibit structure in their absorption spectra. 
The A 2 n 3/2 states of these bands have a short lifetime due to 
spin-orbit coupling between the excited state and one or more 
O = 3/2 repulsive states . 66 A near-cancellation of the overlap 
of bound and continuum wavefunctions at v' = 2 provides a 
comparatively slow predissociation at that level, giving rise to 
narrow rotational absorption features (FWHM ~ 0.04 cm -1 ) 
suitable for LIF excitation. The laser-induced fluorescence ex- 
citation spectrum was measured for this work and is shown in 
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FIG. 2. (a) IO fluorescence excitation spectrum measured in this work and 
(b) the rotationally resolved absorption spectrum of the (2,0) band of the 
A 2 n 3/2 <— X 2 n 3 / 2 transition. 54 


Figure 2 along with the rotationally resolved absorption spec- 
trum of the (2,0) band of the A 2 Fl3/2 X 2 113/2 transition of 
IO reported by Wennberg et al . 24 

Dissociation of IO in the V = 2 level of the A 2 113/2 elec- 
tronic state becomes increasingly competitive with fluores- 
cence as higher J states are excited. The lifetime for the ex- 
cited state of IO is short: ~1 ns for J' = 1.5 decreasing to 
15 ps for J' = 50. 5. 66 Collisional quenching of the V = 2 
rotational states in N2 was characterized by Gravestock 
et al. 6 1 The inferred rate constant of quenching (kq) for the 
bandhead (a mixture of J' = 3.5, 4.5, and 5.5 states) was 
k q = 9.6 ± 6.3 x 10 -11 cm 3 mol -1 s” 1 and increased at higher 
pressures and J' values. Excitation to the J' = 7.5 state of the 
(2,0) band around 445 nm was selected as the basis for detec- 
tion in this work. Using k q from the bandhead, a reasonable 
approximation for the conditions used in this study (J' = 7.5, 


20 Torr) yields a quenching lifetime (r q ) of 16 ns and a fluo- 
rescence lifetime (tf) of 0.5 ns. 

Figure 3(a) shows the excitation and fluorescence scheme 
employed in this work. The spectrally resolved fluorescence 
emission resulting from the excitation of IO at 445 nm 
(20 Torr) was collected using a 1-m focal length multi-pass 
Ebert spectrometer (Jarrell Ash, Boston, MA). The spectrum 
was then corrected for the quantum efficiency of the photo- 
multiplier tube (PMT) at different wavelengths and is shown 
in Figure 3(b). The fluorescence signal corresponds to the vi- 
brational bands of the A 2 Fl3/2 (v' = 2) -> X 2 Fl3/2 (v" = 0-1 1) 
transitions, confirming that vibrational relaxation is a negligi- 
ble channel compared to fluorescence. The relative intensities 
of the vibrational bands agree well with published Franck- 
Condon factors. 68 In the IO instrument described here, a 
530 nm bandpass filter (Semrock FF01-530/23-25, Rochester, 
NY) and 514 nm longpass filter (Semrock LP02-514RE-25, 
Rochester, NY) are used to capture the fluorescence from the 
v' = 2 -* v" = 5 and v" = 6 transitions as indicated in 
Figure 3(b). These two bands constitute 16% of the total flu- 
orescence intensity arising from the v' = 2 state. 

To maximize the specificity of the IO measurement, the 
laser is tuned across the R7.5 line (445.033 nm, “online”) and 
then to an off- resonance wavelength (“offline”) between the 
R6.5 and R7.5 absorption features (445.028 nm). The differ- 
ence between the signals at the two frequencies is used to de- 
termine the IO mixing ratio. The online frequency has a large 
differential cross section (2.92 xlO” 16 cm 2 mol” 1 ) and the 
offline frequency is characterized by a broad, clean baseline 
at the pressure used for detection (20 Torr). The online and 
offline frequencies were selected to minimize the differential 
absorption cross section of NO2, which has structured absorp- 
tion features in this wavelength region and can be present at 
high concentrations in the atmosphere. The differential cross 
section between the online and offline wavelengths is small 
(<3%) for NO2. Additionally, since the fluorescence of NCB 



FIG. 3. The LIF detection of IO. (a) IO is pumped via a ro-vibronic transition of the A 2 n 3 /2 (v' = 2) X 2 n 3 /2 (v" = 0) band near 445 nm. Red shifted 

fluorescence is observed, (b) The peaks correspond to the vibrational bands of the A 2 n 3 /2 (v' = 2) — >■ X 2 n 3 /2 (v" = 0-11) transitions. Vibrational energy 
transfer to v' = 0 or 1 was not observed. The fluorescence from the A 2 n 3 /2 (v' = 2) — y X 2 II 3 /2 (v" = 5-6) transitions is collected for the IO instrument using a 
combination of a bandpass and longpass filter to reject short- and long-wave interferences. 
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TABLE I. Differences between Harvard and Leeds LIF Instruments. 


Instrument feature 

This work 

Whalley et al . 8 

Max laser power at 445 nm 

80 mW 

150 mW 

Detection axis 

32-pass white cell 

Single pass cell 

Online detection strategy 

Tuning across R7.5 line 

Online acquisition at a single frequency 

Pressure of detection 

20 Torr 

150 Torr 

Obs. fluorescence bands 

(2,5) and (2,6) bands 

(2,5) band 

Reference cell 

Yes 

Not during deployment 

Sensitivity 

0.3 pptv in 60 s 

0.3 pptv in 300 s 


occurs over a broad range of frequencies, the use of long- 
pass and band-pass filters significantly reduces its contribu- 
tion to the count rates observed both online and offline. No 
detectable contribution of NCb interference to the IO signal 
was observed during laboratory testing of up to 300 ppbv of 
NO 2 . Interferences from chemical pathways that could pro- 
duce IO within the instrument were also explored and found 
to be negligible (see Sec. SI of the supplementary material 69 ). 

III. INSTRUMENT DESCRIPTION 

The IO instrument was designed as a first generation 
ground-based/mobile sensor that is automated and integrated 
into an electronic architecture compatible with a broad range 
of field deployment methods. This section provides a detailed 
description of the hardware. Numerous considerations had to 
be addressed during the design process: 

Sampling. The optimal design would draw sampled air 
to the detection axis in an essentially “wall-less” manner. IO 
radicals are extremely reactive and are lost efficiently to in- 
ternal surfaces. A reduced pressure detection axis is used in 
which the distance between the inlet and detection region are 
minimized to meet these criteria and minimize losses. 

Laser subsystem. A narrow-bandwidth tunable visible 
laser is required. High pulse repetition rates with average 
power greater than 20 mW are needed to make high signal-to- 
noise measurements of the sub-pptv mixing ratios of IO in the 
atmosphere. The laser subsystem is a nanosecond Nd:YAG 
pumped Titanium: Sapphire laser with a compact footprint. 

Detection axis. The fluorescence detection axis was de- 
signed to allow large throughput while heavily discriminating 
against solar scatter and laser scatter. The detection axis se- 
lected is a 32-pass White cell with optical baffling and spectral 
filtering analogous to the Harvard NOi instrument detection 
axes. 60 

Data acquisition. The ground-based version of the IO 
instrument is expected to work for hours with minimal hu- 
man interaction. Programs for data acquisition and experi- 
mental control were written in LabVIEW (National Instru- 
ments, Austin, TX) with error recovery capabilities. 

Future integration. The IO instrument was designed with 
the capability to be integrated with the existing Harvard CIO 70 
and BrO 71 flight instruments. As a result, the detection duct 
is identical to the halogen instruments used for airborne mea- 
surements by our group. 

The instrument described in this work draws on simi- 
lar fundamental spectroscopy as the LIF instrument deployed 


previously by the University of Leeds. s ' 10,48 Major differ- 
ences between the instruments are emphasized in Table I. The 
focus of the design process for this instrument was to increase 
specificity, time response, and mobility to enable sub-pptv 
measurements of IO with spatial and temporal resolution not 
obtained previously. 

A. Sampling 

IO is sampled using the fluorescence assay by gas expan- 
sion (FAGE) technique, a common technique for sampling 
free radicals in the atmosphere. " The sampling assembly is 
presented in Figure 4. The inlet consists of a flat plate with 
a 0.75 mm pinhole aperture for baric filtering and a 5.08 cm 
x 5.08 cm x 35 cm long duct, which delivers the reduced- 
pressure sample to the detection axis. A scroll vacuum pump 
(Varian Triscroll 600, Palo Alto, CA) is used to draw the sam- 
ple through the inlet and into the detection axis. Flow rates 
of the sample through the system are typically 5 standard 
liters per minute. Laboratory tests of the pinhole inlet were 


FIG. 4. Airflow within the IO instrument. Air enters the duct through a pin- 
hole aperture, which reduces the pressure of the sample to 20 Torr. The sam- 
ple then flows directly to the detection axis where a White cell, photomul- 
tiplier tube, optical filters, and an optical light trap are used to detect the 
fluorescence signal. The sample is exhausted through a vacuum pump. 
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conducted with a modular cavity ringdown detection axis 
placed in front of the pinhole aperture. These experiments en- 
abled a direct measurement of the transmission efficiency of 
IO as it passed through the inlet to the detection axis of 83% 
± 4%. The transmission efficiency of the pinhole aperture was 
higher than that of a molecular beam skimmer of the same di- 
ameter (77% ± 5%), which is another commonly used critical 
orifice in FAGE instruments. 

B. Laser subsystem 

A tunable, narrow-linewidth Ti-Sapphire laser (Photon- 
ics Industries DC150-532/Tu-H, Bohemia, NY) was custom 
built as the excitation source for the IO LIF instrument. A 
schematic of the optical design of the Ti: Sapphire system is 
presented in Figure 5. A Q-switched Nd:YAG laser at 532 nm 
(3.2 W average power at 5 kHz) pumps the Ti:Sapphire cav- 
ity. The wavelength of the laser is selected by tuning the angle 
of a diffraction grating (2000 lines/mm). The fundamental of 
the Ti:Sapphire (~890 nm) is doubled with a lithium tribo- 
rate crystal to produce the 445 nm excitation wavelength used 
in this experiment. The TfiSapphire system produces an aver- 
age power at 445 nm of approximately 80 mW at 5 kHz with 
40 ns pulse widths. The high repetition rate of the laser yields 
relatively low per pulse energies (<1.6 x 1 0 s J), reducing 
the likelihood of saturation in the excitation of IO. The up- 
per limit of the linewidth of the Ti: Sapphire, as determined 
from laboratory measurements of IO and NO 2 absorption 
spectra at 20 Torr, was less than 0.05 cm -1 , consistent with 
the theoretical linewidth of 0.03 cm' 1 . The laser linewidth is 
well matched to the linewidth of the R7.5 absorption feature 
(0.04 cm' 1 ) used for the detection of IO. The Ti:Sapphire and 
Nd:YAG pump laser have been integrated into a single, com- 
pact optical package with a 32 cm x 32 cm x 13 cm volume. 



FIG. 5. Schematic of the pulsed visible Ti:Sapphire laser system. The 
Ti:Sapphire laser is powered by a single frequency-doubled Nd:YAG laser. 
The 890 nm output is then frequency-doubled to 445 nm. The footprint of the 
laser system is 32 cm x 32 cm. 


Several modifications to the laser system have been made 
to improve its stability in the field and to aid in the analy- 
sis of laser performance. First, thermo-electric cooling mod- 
ules (Laird Technologies CP2-31-06-L1, London, England) 
have been installed to improve thermal stability of the laser 
baseplate. Second, custom-built silicon photodetectors were 
implemented to monitor the output power of the Nd:YAG at 
532 nm and the fundamental of the Ti:Sapphire at 890 nm. 
Third, a solenoid-driven laser shutter was incorporated to 
collect laser-free background counts during data acquisition. 
Finally, a piezo-driven rotation stage (Physik Instrumente 
M-035.ps, Germany) was implemented in place of a step- 
per motor for wavelength tuning. The piezo-driven rotation 
stage improves the speed and reliability of wavelength tuning. 
The minimum increment of wavelength tuning with the piezo- 
driven stage is 0.01 cm' 1 . A custom-designed circuit provides 
stable, filtered voltages to the piezo element and monitors the 
position of the piezo with high sensitivity measurements from 
the integrated strain gauge sensor. This strain gauge signal 
serves as a convenient metric for the optical frequency of the 
laser. 

C. Detection axis 

The Ti:Sapphire laser output is directed to the detection 
axis by several turning mirrors on the optical breadboard. A 
multi-pass White cell, depicted in Figure 6(a), is configured 
so that the laser makes 32 passes through the detection vol- 
ume before exiting. The White cell mirrors (Spectrum Thin 
Films, Bohemia, NY) are spaced at 15.5 cm with a reflectiv- 
ity of >99.99% at 445 nm. The individual laser beams are 
expanded as they pass through the detection volume and are 
non-overlapping in space (see Figure 6(a)), reducing the po- 
tential for saturation in the absorption step. The entrance mir- 
ror of the White cell is sealed with an AR-coated glass win- 
dow. Typically 90% of the incident laser power is transmitted 
through the cell. A small purge of nitrogen (~ 1 00 standard 
cubic cm per minute (seem) per mirror) is introduced at the 
entrance and twin mirrors to maintain mirror cleanliness dur- 
ing sampling. The detection axis is baffled to reduce inter- 
ference from laser and solar scatter. Knife-edge baffles line 
the walls and the region between the White cell mirrors and 
the detection volume. The baffles and the chamber walls are 
painted with a black paint selected to minimize fluorescence 
(Zuel, St. Paul, MN), and an optical trap is situated perpen- 
dicular to the propagation of laser light and directly opposite 
the detector to reduce chamber scatter. 

A series of lenses and filters (Figure 3(b)) are used to 
direct fluorescence from the 1 .5 cm x 2 cm x 3 cm volume 
in the center of the detection axis onto the gallium arsenide 
phosphide photocathode of the PMT (Hamamatsu Photonics 
H7421-40, Japan), which operates in photon-counting mode. 
On average, the solid angle collected by this lens system is 
2.5% of Arc . The PMT was selected due to its high quantum 
efficiency (~38%) at the fluorescence collection frequencies 
and somewhat lesser quantum efficiency (25%) at the laser 
excitation wavelengths. The PMT is actively cooled to 0°C 
with thermoelectric coolers to reduce dark counts (typically 
<20 counts per second). 
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FIG. 6. A schematic of the optical systems, (a) A 32-pass White cell is used. 
Two lenses image the fluorescence from the central 1.5 cm of the White cell 
onto the photocathode of the PMT inside a refrigerated housing. Filters reject 
signal at the laser wavelength and provide high transmission at the observed 
fluorescence frequencies. Background noise and scattering are suppressed us- 
ing optical baffles and a trapezoidal light trap across from the detector. Gas 
flow is in the direction perpendicular to the plane of the drawing, (b) The 
laser, optics, and detection axes are all configured on a 56 cm x 60 cm op- 
tical plate. The use of multiple vertical planes minimizes the footprint of the 
detection axis and enables each of the three detection axes (White cell, refer- 
ence cell, CRDS cell) to be aligned independently. 

D. Frequency reference cell 

Changes in temperature, pressure, and humidity of air 
within the Ti:Sapphire cavity can result in shifts in the fre- 
quency of the laser on the order of the linewidth of the IO 
absorption features at 20 Torr. Rather than attempting to mon- 
itor the frequency of the laser directly, a reference cell in 
which IO is continuously generated is used to provide a sta- 
ble metric of optical frequency. The frequency reference cell 
is a baffled, pressure-sealed, single-pass cell. The cell is filled 
to 20 Torr with a mixture of 20% CF3I in ultra zero air in 
order to reproduce the line shapes anticipated in the White 
cell detection axis, which is operated at the same pressure. 
Current (~1A) is applied through a Nichrome filament to 
generate IO from the precursor molecules. While the chem- 
ical mechanism is not well defined, it is likely that I atoms 
are generated from CF3I and react with oxygen to produce 
IO. Approximately 12-14 mW of the Ti:Sapphire laser light 


is directed into the frequency reference cell using a neutral 
density filter (Thorlabs NE01BA, Newton, NJ). A photomul- 
tiplier tube (Electro-Mechanical Research 541N-09, Sarasota, 
FL) in gated photon-counting mode is situated perpendicular 
to the axis of laser propagation to detect the laser-induced flu- 
orescence signal at X > 500 nm. A longpass interference filter 
(Semrock LP02-514RE-25, Rochester, NY) is used to reject 
scatter from the laser and other short-wave interferences. The 
concentration of IO in the cell is not quantified because the 
cell serves only as a frequency reference and does not play 
any further role in the calibration of the detection axis. 

The leak rate of ambient air (<1 Torr per day) into the 
reference cell limits the amount of time it can operate prior to 
refilling (~10 days). As the pressure in the cell increases, the 
spectroscopic features broaden. Since the width of the refer- 
ence cell peak is incorporated into an algorithm used to an- 
alyze data from the detection axis, the reference cell is re- 
filled when pressures exceeds 28 Torr. The reference cell is 
equipped with solenoid valves to enable fast evacuation and 
recharging of the gas mixtures, limiting the system downtime 
to refill the reference cell to less than 10 min. 


E. Optical design 

A schematic of the optical design for the IO instrument 
is presented in Figure 6(b). The multilevel system is based 
on a standard optical breadboard with a 56 cm x 60 cm foot- 
print. A series of lenses and turning mirrors condition the laser 
beam and direct it to the White cell, frequency reference cell, 
and a modular cavity ringdown detection axis used for calibra- 
tions. Two silicon photodiodes are used to monitor the laser 
power at 445 nm prior to and upon exit from the White cell 
as a metric for the alignment of the detection axis and to nor- 
malize the fluorescence signal. The compact design enables 
free-space coupling of the laser-beam, which is essential for 
maintaining a high beam quality for the multi-pass White cell 
and cavity ringdown detection systems. 

F. Data acquisition 

The data acquisition system for the instrument combines 
software written in LabVIEW, two 16-bit A/D cards (National 
Instruments PCI-6229, Austin, TX), and a custom-designed 
instrumental distribution panel. This data acquisition system 
is central to the operation of the instrument, controlling the 
commands sent by the software to hardware systems (e.g., 
sending voltages to the grating rotation stage, flow controllers, 
or laser shutter) and the acquisition of data coming from in- 
strumental hardware (e.g., reading signals from PMTs, ther- 
mistors, photodiodes). Telemetry and signal data are written 
continuously at 1 Hz and 3 Hz respectively to a binary file. 

The fluorescence signals from the White cell and ref- 
erence cell are acquired using gated photon counting. The 
pulsed outputs from the PMTs are shaped with a pulse ampli- 
fier/discriminator and buffered with circuitry implemented on 
the instrumental distribution panel. The fluorescence signal 
must be collected during the laser pulse (30-40 ns) due to the 
short fluorescence lifetime of IO (<1 ns). The gated counting 
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FIG. 7. Diagram of the gating scheme employed for photon counting. The 
trigger used to initiate gating occurs at time = 0 s. The measured fluores- 
cence and laser profiles are shown. A fluorescence gate 300 ns long is used to 
capture the signal. A background gate 30 /xs long detects continuous sources 
of background that occur after the fluorescence and short-lived interferences 
have subsided. 
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the IO spectrum. A peak-fitting algorithm is applied to the 
reference cell spectrum, and a logic-based process is used to 
determine the position of the R7.5 rotational line used for the 
measurement. The laser is tuned to the offline position prior 
to the peak, where background counts are recorded for 15 s. 
Online data are then acquired by tuning over approximately 
the top 50% of the R7.5 line in approximately 15 s. This on- 
line scan is repeated 15 times. After each online scan, a peak 
fitting analysis is performed to determine the peak center and 
adjust the scan range. Data are then acquired at the offline po- 
sition to collect background counts for an additional 15 s. The 
online and offline acquisition sequence is repeated 15 times 
before reacquiring the IO spectrum. Status flags are assigned 
to each of these steps and recorded in the data stream to facil- 
itate data analysis. In practice, this method of data acquisition 
is very reliable since the peak-fitting procedures can correct 
for any drift in laser frequency. 


reduces the influence of solar scatter and other long-lived flu- 
orescence interferences. A schematic of the gating scheme is 
presented in Figure 7. To generate the gate, the pulsed out- 
put from the Q-switch of the Nd:YAG is used as a trigger. 
A counter on the LabVIEW A/D card generates a 300 ns long 
gate delayed 1 /xs from the rising edge of the Q-switch trigger. 
The “pause-trigger” functionality within LabVIEW is used to 
count the pulses from the PMT that arrive within the gate. For 
the White cell only, a second gate is used to collect the back- 
ground signal that arises from solar scatter and dark counts. 
The background gate is 30 /xs long and is delayed 2.4 /is 
relative to the rising edge of the Q-switch trigger. The gated 
counts are integrated over 330 ms (1650 laser pulses) before 
being recorded. 

The software interface for the system is designed to op- 
erate autonomously for extended periods of time with no in- 
tervention. During normal operation, the system repeats a se- 
quence of events to collect and record data. Figure 8 shows 
an example of the raw data collected during a typical acqui- 
sition sequence. First, the laser is physically shuttered so that 
dark counts can be recorded. The laser is then scanned from 
445.00 nm to 445.05 nm across several rotational features in 



FIG. 8. The data acquisition sequence for the IO instrument. Off-resonance 
(cyan dots), and shuttered (black dots) measurements are collected. The laser 
is then tuned across several rotational lines in the IO spectrum (red line) and 
then is repeatedly tuned across the R7.5 rotational line (green). The R7.5 line 
is indicated within the black dashed rectangle in the IO spectrum and in green 
representing online data acquisition. 


G. Deployment and operation 

A mobile, weatherproof platform was developed for field 
deployment of the IO instrument. The instrument was de- 
signed to be compatible with a broad range of deployment 
conditions, particularly marine and polar environments where 
halogen chemistry is most active. An instrument deployed to 
these environments must be capable of withstanding tempera- 
ture fluctuations, wind, precipitation, and corrosion. In its de- 
ployable package, the IO instrument is divided into two parts: 
(1) the sensor unit, which contains the inlet, laser, detection 
optics, and associated data acquisition electronics; and (2) 
the auxiliary unit, which contains the vacuum pump, ground 
fault circuit interrupters, uninterruptible power supply (UPS), 
and a chiller or water circulating pump. The elements in each 
unit are housed in a 1.3 m x 0.7 m x 0.6 m polypropylene 
copolymer case (Pelican Products Inc. 0550, Torrance, CA) 
customized with weather-tight plumbing and electrical con- 
nectors as well as the addition of brackets and vibration iso- 
lators to improve the stability of the optical alignment. The 
cases are fitted with chassis and large, low-pressure tires de- 
signed for movement over rugged terrain. The units can be 
configured not only with handles for movement by hand, but 
also have a ball hitch to enable towing with an all-terrain ve- 
hicle over longer distances. 

Two different cooling strategies were designed for tem- 
perature control within the sensor unit. In both schemes, the 
sensor unit is cooled by a liquid loop that flows in parallel 
through two fan-circulated heat exchangers and the laser base- 
plate and diodes. Temperature sensors (E5CK-AA 1-500, Om- 
ron, Kyoto, Japan) control the fans to prevent extreme fluctu- 
ations in operating temperature. In the first cooling scheme, 
a recirculating chiller (Thermo Neslab Inc. Coolflow CFT-75, 
Newington, NH) supplies the liquid loop with water of an ap- 
proximately constant temperature (18 °C). This strategy was 
used in the field as well as in laboratory preparation of the 
instrument. The alternative strategy harnesses passive cooling 
by the ocean. A 30 cm x 20 cm aluminum heat exchanger 
is submerged in the ocean and a pump is used to circulate 
water through the heat exchanger and sensor box in a closed 
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loop. This approach has the advantage of requiring less energy 
and eliminating the need for refrigerants in sensitive environ- 
ments. However, this strategy requires that the water pump be 
adequately sized to account for the vertical distance between 
the heat exchanger and the instrument, which depends on the 
operating environment. 

The power requirements of the instrument are significant, 
~ 4 kW peak power, with the majority of the power consumed 
by the scroll pump and chiller. The instrument is powered us- 
ing 230V AC. In field environments without adequate local 
power supplies, a 10 kW diesel generator (Kubota Corpora- 
tion GL 1 1000TM, Osaka, Japan) accompanies the instrument 
and can be situated approximately 100 feet away from the sen- 
sor unit. The generator is sized to ensure that the entire system 
can operate for more than 12 h without refueling. 

IV. DATA ANALYSIS: CALCULATION OF THE MIXING 
RATIO OF 10 

Post-acquisition processing of the data is required to con- 
vert the observed signals to atmospheric mixing ratios of IO. 
This section will describe the data processing steps includ- 
ing: normalizing the IO data, performing a nonlinear fit of the 
spectral feature, analysis and internal consistency tests to re- 
move potential interferences, and tying the observed signal to 
IO concentrations by establishing an absolute calibration. 


the shoulders, fitting the line to a Gaussian function is a rea- 
sonable and computationally efficient approximation. 

The first step in the fitting procedure is to perform a non- 
linear least-squares fit of the reference cell data, which have 
high signal to noise ratios. The fits of the reference cell peaks 
show high R 2 values (>0.98) and no significant structure in 
the residuals, confirming that the Gaussian model is a good 
approximation for the observed line. The line center deter- 
mined in the reference cell is then used to constrain the line 
center for the Gaussian fit of the White cell data, which typ- 
ically have lower signal to noise. The quality of the fit is as- 
sessed using the goodness of fit parameter and the linewidth. 
These constraints are used to identify fitting errors that may 
result from the successful fitting of optical noise. The height 
of the peak derived from the fit is the proxy for the IO fluo- 
rescence signal, Sio, in counts V~* s -L 

Sio can be related to the mixing ratio of IO in the atmo- 
sphere [IO] atm with the following equation: 


U 0] atm 


Sio 


T f C 


LIF 


(3) 


where TE is the transmission efficiency (0.83 ± 0.04) deter- 
mined from laboratory tests of the inlet, and Clif is an ex- 
perimentally determined calibration factor (counts s -1 V" 1 
pptv -1 ), which will be described in Sec. IV C. 


A. Calculating background-free, 
power-normalized signal 

Continuous spectroscopic interferences and the impact 
of fluctuations in laser power must be removed from the ob- 
served fluorescence signal. First, the count rate measured dur- 
ing the background gate is subtracted from the count rate dur- 
ing the fluorescence gate. The observed signals from White 
cell and reference cell are then divided by the measured laser 
power (in Volts) to account for the linear dependence of flu- 
orescence and other laser-induced interferences (e.g., laser 
scatter) on laser power. Finally, the mean of the offline data 
segments collected before and after each online peak scan is 
subtracted from the online peak scans. 

B. Peak fitting 

After the fluorescence signal of IO is isolated from the 
background, it is quantified using a peak fitting procedure. 
Since the laser is tuned across the R7.5 rotational line in the 
ro-vibronic spectrum, the online signal reflects the lineshape. 
For the specific conditions of this measurement, the lineshape 
is dictated by a combination of Lorentzian terms, i.e., natural 
broadening (HWHM = 2.33 x i 0 4 nm) and pressure broad- 
ening (HWHM = 4.25 x 1 0 6 nm), and Gaussian terms, i.e., 
Doppler broadening (FWHM = 4.586 x 1 0 4 nm) and laser 
broadening (FWHM = 6.62 x 10 -4 nm). While a Voigt pro- 
file is the best model for a lineshape that is a convolution 
of Gaussian and Lorentzian terms, the fitting procedure for 
a Voigt profile is computationally expensive. Since the data 
acquisition system narrows the range of wavelengths sampled 
to emphasize the center of the spectral line at the expense of 


C. Calibration 

The fluorescence signal is related to a mixing ratio of 
IO by performing a standard addition experiment with cali- 
brated amounts of IO. IO is generated in the laboratory by 
titrating known quantities of a certified standard of 9 nmol 
mol -1 CF3I in N2 (Scott-Marrin, Inc., Riverside, CA) with an 
excess of atomic oxygen. The CF3I is diluted with a bulk flow 
of N2 to maintain the overall flow rate and is introduced into 
the 32-pass LIF detection axis through a calibration assembly 
that is mounted on top of the instrumental inlet. Atomic oxy- 
gen is generated by passing a 0.6% mixture of O2 in helium 
through a microwave discharge (Scintillonics Inc. HV15A, 
Fort Collins, CO). A custom designed quartz loop injector 
(James Glass, Hanover, MA) is used to introduce atomic oxy- 
gen to the center of the flow approximately 4 cm downstream 
of the pinhole inlet where it mixes with the CF3I/N2. IO is 
produced from the following reaction with a known yield of 
83% ± 11%: 73 

CF 3 I +0^-10 + CF3. (4) 

CF3I is added by standard addition with a constant flow 
rate of O2. The experiment is then repeated over a range of O2 
flow rates to determine the conditions under which the pro- 
duction of IO is limited by available CF3I rather than atomic 
oxygen. Figure 9(a) shows the fluorescence signal as a func- 
tion of O2 flow rate with a constant amount of CF3I added. 
The plateau in fluorescence signal at higher flow rates of O2 
that is observed suggests that O is present in excess under 
those conditions. The results of a standard addition experi- 
ment conducted with 2 different O2 flow rates are depicted in 
Figure 9(b). 
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FIG. 9. Results from a typical chemical titration calibration. This experiment was performed at the Shoals Marine Laboratory on August 18th, 2011. (a) A 
standard addition experiment is performed at a variety of O 2 flow rates, shown here for a constant 400 seem flow of 9 ppbv CF 3 I in N 2 . An excess of atomic 
oxygen indicated by the plateau in signal change at higher O 2 flow rates, (b) A standard addition of IO is performed with a constant concentration of added O 2 . 
The slope of the linear fit of the standard addition data is the calibration factor. 


The concentration of IO measured at the detection axis 
is determined from the concentration of added CF3I, and the 
yield of Reaction 4 (83% ± 1 1%). A kinetic model is used to 
account for the chemical losses that occur between the time 
of production of IO and its detection. The dominant chemical 
loss pathways and associated rate constants 74 in units of cm 3 
mol -1 s -1 are: 

70 + 70 OIO + I 

k = .44 x (1.5 x IO -11 ) x exp(500/T), (5) 

70 + 70 27 + 0 2 

k = .56 x (1.5 x KT 11 ) x exp (500/ T) , (6) 

70 + 0-^7+02 k = 1.2 x 1(T 10 . (7) 

Since the calibration is performed in excess atomic oxygen, 
the losses due to Reaction 7 are significant. Therefore, a di- 
rect measurement of the concentration of atomic oxygen is 
required to quantify chemical losses in the kinetic model. 

The concentration of atomic oxygen is measured by in- 
corporating an absorption axis directly downstream of the po- 
sition where atomic oxygen atoms are introduced into the 
flow. The light source is an atomic oxygen lamp, fabricated in 
our laboratory with a KMnC>4 source. The absorption of the 
oxygen line ( 3 S 1 •<— 3 Po) at 130.6 nm is monitored with a 0.2 
meter vacuum monochromator (Acton Research Corporation 
VM-502, Acton, MA). Using the Beer-Lambert law, the con- 
centration of atomic oxygen is calculated from the measured 
transmittance, the absorption cross-section (1.64 x 10 -13 
± 0.36 x 10 -13 cm 2 mol -1 5 ), and the pathlength (5.08 cm). 
The yield of atomic oxygen relative to the molecular oxy- 
gen added is low (<5%), which is likely due to a combina- 
tion of inefficient production of O by the discharge and large 
losses within the quartz loop injector. The concentration of O 
is incorporated into the chemical loss model to correct for the 


losses of IO that occur within the duct. A 22% uncertainty in 
calculated concentrations of O is derived from uncertainties in 
the absorption cross section (20%) and absorption measure- 
ment (11%). Above a threshold level of 25 ppmv of added 
O2, the yield of IO produced did not show significant depen- 
dence on additional O2 for the conditions of the calibrations 
performed in this work. Only calibration data performed at 
O2 mixing ratios greater than 35 ppmv and less than 60 ppmv 
were used in sensitivity calculations to ensure that small vari- 
ations in the concentration of atomic oxygen produced would 
not affect the result. 

The sensitivity of the LIF detection is determined from 
the slope of a least squares regression of the fluorescence sig- 
nal as a function of the IO mixing ratio introduced into the 
system (e.g.. Figure 9(b)). The observed linearity in LIF sig- 
nal observed as a function of IO mixing ratio implies that the 
absorption step is not saturated and the PMT is operating in a 
linear range over atmospherically relevant IO concentrations. 
The sensitivity of the LIF detection is reported as a calibration 
factor, Clif, in units of counts s _1 (pptv IO) -1 V 1 at 20 Torr. 
The standard addition experiments in Figure 9 were per- 
formed during a field deployment of the IO instrument. Two 
different concentrations of added O2 yielded Clif values from 
the fitting procedure of 25.0 ± 2.0 and 24.8 ±2.1 counts s -1 
(pptvIO)- 1 V" 1 . 

The instrument has also been designed to incorporate a 
secondary method of calibration by cavity ringdown spec- 
troscopy (CRDS), which has orthogonal sources of uncer- 
tainty relative to chemical titration. A modular CRDS axis 
is configured downstream of the LIF detection axis so that 
calibration can be tied to the absolute absorption cross 
section.'’ 4 - 76 This method is less appropriate for field cali- 
brations. Given the pathlength of the cell (22 cm) and mir- 
ror reflectivity (currently 99.95% at 445 nm), concentrations 
of IO above 8 x 10 8 molecules of IO per cm 3 (~1 ppbv 
at 20 Torr) are required for CRDS detection. The introduc- 
tion of high concentrations of IO precursors can contaminate 
the surfaces within the instrument and reduce the reflectivity 
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of the mirrors. While the CRDS calibration was not im- 
plemented in the field, calibrations were performed previ- 
ously with higher reflectivity (>99.99%) mirrors (Los Gatos 
Research — no longer available, Mountain View, CA) in the 
laboratory. Calibration factors determined with CRDS and 
LIF in the laboratory agree to within less than 10%, well 
within the uncertainty of the individual methods. 

D. Uncertainty 

Given the relatively low concentrations of IO anticipated 
in the atmosphere, a thorough analysis of the uncertainty of 
the IO measurement is required. To quantify the reproducibil- 
ity of the peak fitting technique used to determine the peak 
height, a statistical bootstrap analysis of the fitting procedure 
was conducted to estimate the variance. This analysis is de- 
tailed in Sec. S2 of the supplementary material. 69 The calcula- 
tion of the yield of IO from CF3I (11%), the laser power mea- 
surement (3%), and the fitting procedure (8%) all contribute to 
an overall 2 a measurement uncertainty of 14%. The end-to- 
end estimate of measurement uncertainty for the instrument is 
21% (2er), resulting from the 14% and 16% uncertainties in 
the determination of Clif and TE, respectively. 

The uncertainty at low concentrations of IO is dominated 
by the peak fitting since the other sources of uncertainty are 
relative and scale with concentration. The bootstrap analysis 
revealed that the 1 a uncertainty associated with the fitting 
procedure is 0.2 pptv for each individual peak fit. Four peak 
fits are averaged for each reported 1-min data point, so the la 
uncertainty, which scales with i, where N = 4, is 0.1 pptv 
for the minute-averaged measurement. The limit of detection, 
defined as the threshold concentration of IO where the signal 
to noise ratio equals 3, is 0.3 pptv in 1 min. 

It should be noted that the excitation rate of IO de- 
pends on temperature. The use of a narrow-linewidth exci- 
tation source in this experiment enables the excitation of elec- 
trons from a single rotational level (v" = 0, J" = 7.5) in 
the X 2 Il3/2 state of IO. The normalized population density 
of the rotational state will vary as a function of temperature. 
The calibrations and measurements presented in this work oc- 
curred within a relatively narrow temperature range since the 
detection axis temperature was tightly regulated. However, 
corrections to the calibration using ambient temperature ob- 
servations would be required for airborne measurements. The 
relative population densities can be calculated directly using 
spectroscopic constants and measurements of temperature. 59 

V. FIELD VALIDATION 
A. Site description 

The prototype instrument was deployed to Shoals Marine 
Lab (SML) on Appledore Island, ME, a marine science fa- 
cility affiliated with Cornell University and the University of 
New Hampshire. Appledore Island (42°58' N, 70°37' W), ap- 
proximately 10 km from Portsmouth, NH, is the largest island 
within an archipelago known as the Isles of Shoals. This site 
was chosen for the validation of the Harvard IO instrument 
because Stutz et al. 40 measured IO by LP and MAX-DOAS 



FIG. 10. A map of the Isles of Shoals. The light path of the long path DOAS 
instrument deployed in 2004 is indicated spanning from Appledore Island 
to White Island. LIF measurement sites are indicated (yellow stars). From 
August 17th to August 25th, 2011 measurements were made at site 1. From 
August 31st to September 5th, measurements were made at site 2. Kelp stands 
had been observed near Siren’s Cove, Devil’s Dance Floor, and Norwegian 
Cove in previous years. (Map adapted with permission from W. E. Beamis, 
Shoals Marine Laboratory.) 

there previously as a part of the Chemistry of Halogens at the 
Isles of Shoals (CHAiOS) field campaign in 2004. The light 
path of the LP-DOAS instrument extended from Appledore 
Island to neighboring White Island (Figure 10) and ranged 
from 15 to 35 m above the ocean surface. Reactive iodine 
species (IO and OIO) were measured by the LP-DOAS in- 
strument on 19 of 23 days, with a maximum mixing ratio of 
4 pptv for IO. 

The Harvard IO instrument was deployed for validation 
from August 17th to September 5th, 2011. Data were not col- 
lected on August 20 due to heavy rains or during the period 
from August 26th to August 30th due to the onset of Hurri- 
cane Irene. Observations were made from two different mea- 
surement sites on opposite sides of the island, as depicted in 
Figure 10. From August 17th to August 25th, observations 
were made from the western shoreline of the island near Sand 
Piper Beach (site 1). This relatively sheltered site was se- 
lected due to a combination of prevailing onshore winds and 
the accessibility of terrain near the high-tide line. Observa- 
tions from August 31st to September 3rd were carried out on 
the eastern shoreline of the island, near Broad Cove (site 2). 
Broad Cove has sheer walls and is significantly more exposed 
than Sandpiper Beach. 

Laminariales (kelp) considered largely responsible for 
the tidal signatures observed in emissions at Mace Head, 
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Ireland and Roscoff, France, were not observed to be actively 
growing at either site. 10, 13,28 However, abundant Saccharina 
latissima and limited quantities of Laminaria digitata washed 
ashore in the days after Hurricane Irene near Broad Cove. In- 
dividuals familiar with the ecology of the island indicated that 
kelp stands grow most frequently on the northern edge of the 
island near Siren’s Cove as well as in the area surrounding 
Norwegian Cove. 77 These northerly sites were inaccessible 
for sampling due to densely wooded terrain. 


B. Field performance 

The sensitivity of the measurement is directly impacted 
by changes in laser power and fluctuations in the alignment 
of the White cell. Laser power decreased by ~40% early in 
the deployment. Humidity changes within the sensor box and 
resulting condensation onto optics within the laser was likely 
the cause of the decreased power. To prevent condensation, 
the sensor unit, and the laser specifically, were aggressively 
purged with nitrogen and packets of a silica-based drying 
agent were used to keep the relative humidity below 25%. 
Dry operating conditions reversed laser power losses and in- 
creased the laser stability for the remainder of the deploy- 
ment. The practice of normalizing signal by observed laser 
power mitigates the impact of fluctuations in power on the 
measurement. During the field deployment, the optical sys- 
tem was remarkably stable. No realignment of the White cell 
was required after leaving the laboratory. There were no rapid 
changes in transmission and the overall change in alignment, 
even after transporting the instrument, was less than 10% 
throughout the deployment. The average transmission of the 
laser through the White cell during the deployment was 87%. 
Despite very high transmission through the cell, elevated lev- 
els of background counts were observed during the fluores- 
cence gate. In order to avoid nonlinearity in photon counting 
by the PMT, only a fraction of the total laser light was cou- 
pled into the cell, typically on the order of 20-25 mW. Field 
calibrations were performed to assess the overall sensitivity 
of the instrument. The field calibrations agreed to within 5% 
with calibrations performed in the laboratory. 


C. Results 

The observed mixing ratios of IO during the deployment 
to Shoals Marine Laboratory are presented in Figure 11. It 
should be noted that on DOY 23 1 , replacing the anodized alu- 
minum inlet with a version optically blackened with a copper 
oxide coating (Ebonol-C) significantly reduced the solar scat- 
ter detected by the PMT. While the inlet does not account for 
all scatter within the data, it significantly reduced background 
optical counts and reduced errors in the peak-fitting algorithm 
used to determine the IO mixing ratio. 

Minute-averaged mixing ratios ranged from below the in- 
strumental detection limit (0.3 pptv) up to 10 pptv. The max- 
imum IO signal corresponded to daytime low tide events, 
which is consistent with observations of IO at coastal sites 
in Europe. 7_I " ’ 4,4S This tidal dependence is attributed to the 
exposure of macroalgae to oxidative stress, resulting in the 
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FIG. 1 1 . Timeseries of observed IO mixing ratios, (a) Timeseries of IO mix- 
ing ratio (solid blue circles) and tidal height (green line) observed at the 
Shoals Marine Lab during late August and early September 2011. (b) IO 
(solid blue line) with 2<j uncertainty (light blue envelope) observed at site 2 
on September 1st (DOY 245). 


release of molecular iodine and photo-labile iodocarbons. s 
Figure 12(a) shows the IO signal as a function of time of day. 
A diurnal trend is apparent in the IO observations with the 
lowest mixing ratios observed at night and generally higher 
mixing ratios with significantly higher variability during the 
day. IO was measured during 4 complete nights in the ob- 
servation period. The majority of measurements were at or 
below the detection limit of the instrument, but IO was mea- 
sured above the detection limit on DOY 233/234 with a 
mixing ratio of 1 .4 pptv. These observations took place im- 
mediately following Hurricane Irene, when there was gen- 
erally more scatter in the observations. The positive detec- 
tion may have arisen from large quantities of macroalgae 
that had been washed ashore above the high tide line near 
the instrument. The only known pathway for IO produc- 
tion in the absence of photolysis is the reaction of L with 
NO3. 79,80 Figure 12(b) shows the dependence of IO on tidal 
exposure. Tidal exposure was calculated by taking the in- 
verse of tidal height and linearly mapping it onto a scale 
from 0 to 1, where 1 represents the lowest tide observed 
during the observation period. While maximal mixing ra- 
tios of IO were observed during low tide events (Figure 11), 
the overall correlation between IO and tidal height is low (R 
= 0 . 01 ). 
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FIG. 1 2. Dependence of IO on time of day and tidal height, (a) Binned hourly 
averaged IO signal as a function of time of day (blue circles) shows a clear 
diurnal dependence, (b) IO signal as a function tidal height shows no correla- 
tion (R = 0.01). (c) IO is plotted as a function of the product of tidal exposure 
and solar output (R = 0.65). Tidal exposure and solar output were mapped 
onto a scale from 0 to 1 such that a value of 1 on the plot represents maximal 
solar output and the lowest tides observed during the observation period. 


D. Discussion 

1. Comparison with other in situ observations 

Maximum mixing ratios of IO observed at Shoals Ma- 
rine Laboratory (~10 pptv) were not as large as those ob- 
served by in situ measurements in European coastal sites 
(~30 pptv) with large laminaria populations, 10,48 but the ob- 
servations show similar tidal dependence. A major difference 
between European measurements and those in New England 
is the relationship between the tides and solar output. In Mace 
Head and Roscoff, the lowest spring tides occurred close to 
solar noon. At the Isles of Shoals, the low spring tides oc- 
curred at dawn and dusk. The observations of IO in coastal 
regions are dependent on 2 major factors: (1) the exposure 
of macroalgae to oxidative stress and (2) the generation of IO 
from iodine-containing precursors, which is governed primar- 
ily by photolysis. Accordingly, midday spring tides that ex- 
pose macroalgae that efficiently emit iodine provide optimal 
conditions for the production of IO. At the Isles of Shoals, 
photolysis rates were low during the most extreme low tide 
events, reducing the likelihood of observing high concentra- 
tions of reactive iodine based on currently understood mech- 
anisms. The product of tidal exposure (described in Sec. V C) 
and solar output provides a useful metric for the constraint 
on the production of IO. The photolysis rate constant for CI 2 , 
as calculated from the National Center for Atmospheric Re- 
search Tropospheric Ultraviolet and Visible Radiation (TUV) 
model, 81 is used as the metric of solar output and is normal- 


ized to a scale ranging from 0 to 1 . This metric of solar ex- 
posure is only qualitative because there are differences in the 
photolysis frequencies of CL and I 2 , and the model does not 
account for deviation from clear-sky conditions. The correla- 
tion (R = 0.65) between IO and the product of solar output 
and tidal exposure is significantly higher than for tidal height 
alone. 

Differences in local ecology may also account for the 
lower mixing ratios and the low correlation with tidal height 
reported here relative to the European sites. At both Mace 
Head and Roscoff, there are large laminaria populations in 
the intertidal range, including Laminaria digitata, which is a 
particularly potent emitter of iodine when it is exposed at the 
lowest tides. In contrast, laminaria populations are declining 
in the Gulf of Maine with losses of as much as 95% seen in 
some kelp beds over a 20 year period. 82 No kelp stands were 
observed growing in the intertidal range of Appledore Island 
during the sampling period. Testing of the macroalgae sam- 
ples collected at the site (see Sec. S3 of the supplementary 
material 69 ) revealed that the strongest emissions from locally 
growing species came from Fucus vesiculosus, which was 
abundant in the intertidal range. If Fucus vesiculosus was re- 
sponsible for the emissions, it provides a potential explanation 
the low correlation of IO emissions with tidal height, since it 
is fully exposed during larger portions of the tidal range and 
emits I 2 for a longer duration after it is initially exposed rel- 
ative to Laminaria digitata. 

2. Comparison with LP DOAS observations 
at Appledore Island 

The observations reported here, specifically the magni- 
tudes of IO measured, agree well with the LP-DOAS obser- 
vations at the Isles of Shoals. 411 One difference between these 
observations and the DOAS measurements is that the latter 
indicated no tidal signature. This discrepancy may be due 
to the combination of inherent differences in spatial cover- 
age and the generally weak tidal dependence described previ- 
ously. Only a small fraction of the LP-DOAS light path is in 
the intertidal range. As a result, the LP-DOAS measurement 
may be more sensitive to IO produced from the photolysis of 
longer-lived precursors transported from the intertidal regions 
or from emissions over the open ocean, which are not well un- 
derstood. While the lack of coordinated measurements in this 
study limits the extent of the analysis of this dataset, the mea- 
surements at SML validate the field readiness of the Harvard 
IO instrument and provide a first in situ point measurement of 
IO in North America. 

VI. CONCLUSIONS AND FUTURE WORK 

There are several requirements of an IO instrument suit- 
able for improving scientific understanding. First, the instru- 
ment must be highly sensitive since mixing ratios of IO are 
frequently below 5 pptv in the boundary layer. To achieve 
high precision and accuracy, the instrument must be free of 
interferences and be well calibrated for potential losses of IO 
given its reactivity. Second, the instrument must be capable of 
probing changes in the concentration of IO in space and time. 
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Mobility and fast time response are essential for measuring 
concentrations of IO in a range of environments and altitudes. 

We have built and successfully deployed an instrument to 
measure IO in situ by laser-induced fluorescence that meets 
these criteria. Data from a field deployment to the Shoals Ma- 
rine Laboratory demonstrated the sub-pptv sensitivity of the 
instrument under field conditions. By tuning across the spec- 
troscopic line and implementing a reference cell, the potential 
for interference is reduced and fast time response is retained. 
The multipass detection scheme yielded a 0.3 pptv limit of de- 
tection (3a) in 1 min. Further, we have demonstrated that the 
instrument, which has a compact, mobile, and weatherproof 
design, is adequately stable for deployment on a range of plat- 
forms, including aircraft. While the Harvard 10 instrument is, 
to our knowledge, the most sensitive implementation of an 
in situ measurement reported, several improvements could 
further increase the sensitivity of the instrument and reduce 
the detection limit. Reductions in chamber scatter achieved 
with improved optical coatings within the White cell will en- 
able a greater fraction of the laser power to be used for detec- 
tion yielding increases in sensitivity. 

The IO instrument has been designed to be fundamen- 
tally compatible with aircraft based deployment for broader 
tropospheric and stratospheric coverage. While the instru- 
ment could be deployed on its own, a more potent experi- 
ment would be to combine the IO system with the existing 
BrO, CIO, and NO 2 aircraft instruments built by our labora- 
tory. Integrating the IO instrument with the BrO and CIO de- 
tection axes would provide significant power and weight sav- 
ings since the measurements would share several components 
(i.e., inlet and vacuum pump). The merging of these instru- 
ments would create a payload uniquely capable of probing 
the halogen budget throughout the atmosphere and resolving 
open questions about the role of IO throughout the free tropo- 
sphere and lower stratosphere. 
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